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We report a preliminary value for the zero 
magnetic field Na ^S(f = 1, m = - 1) + 
Na ^S(f= 1, m= — 1) scattering length, 
a\-\. This parameter describes the low- 
energy elastic two-body processes in a di- 
lute gas of composite bosons and deter- 
mines, to a large extent, the macroscopic 
wavefunction of a Bose condensate in a 
trap. Our scattering length is obtained from 
photoassociative spectroscopy with sam- 
ples of uncondensed atoms. The tempera- 
ture of the atoms is sufficiently low that 
contributions from the three lowest partial 
waves dominate the spectrum. The ob- 
served lineshapes for the purely long-range 
07 molecular state enable us to establish 



key features of the ground state scattering 
wavefunction. The fortuitous occurrence 
of a p-wave node near the deepest point 
{Re = 72 ao) of the Og potential curve is 
instrumental in determining a\-\ = (52 ± 5) 
uo and 02,2 = (85 ± 3) oo, where the latter 
is for a collision of two Na ^S(f= 2, m = 2) 
atoms. 
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1. Introduction 



Last year two groups reported the observation of 
Bose-Einstein Condensation (BEC) in dilute gasses of 
ultra-cold ^'Rb and ^^Na [1,2], and another reported 
evidence for reaching the quantum degenerate regime in 
'Li [3] but without observing BEC [4]. The observation 
of BEC in a weakly-interacting gas opens up a whole 
range of possibilities, from fundamental studies of the 
coherent atomic samples produced, to the construction 
of the atom-analog of a laser. Theoretical descriptions 
of the weakly interacting Bose condensate are only now 
being developed and experimental techniques to probe 
the condensate are just beginning to be explored. 
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One of the fundamental parameters required to under- 
stand the approach to BEC and the properties of the 
condensate is the 5 -wave scattering length. This scatter- 
ing length determines the low energy elastic scattering 
rate and thus the evaporative cooling rate as well as the 
nonlinear coupling parameter in the Gross-Pitaevski 
equation [5] for the condensate wavefunction. It is not 
necessary to produce a condensate to measure the s- 
wave scattering length: temperatures in a magneto-optic 
trap (IVIOT) are sufficiently low (~ 1 mK) to limit scat- 
tering to a few partial waves and thus permit a determi- 
nation of the 5 -wave scattering length. 

We probe the scattering wavefunction using the tech- 
nique of photoassociation spectroscopy [6-10]. Two Na 
atoms colliding along the ground state 3^S + 3^S poten- 
tial can absorb a photon to produce a bound molecule. 
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in our case to vibrational levels with energy near the 
3^S + 3^P3/2 asymptote. We detect the formation of 
molecules by sending in a second photon which excites 
the molecule to an autoionizing state, thereby producing 
an easily detected NaJ ion. The relative intensities of the 
molecular photoassociation lines carry information 
about the ground state wavefunction. In particular, we 
find that two specific rovibrational lines that arise from 
p -wave scattering are significantly weaker than the cor- 
responding lines for other nearby vibrational levels. This 
indicates that the former rovibrational state is centered at 
an internuclear separation near a node in the /j-wave 
ground state wavefunction. With the location of this 
node established, the intensities and lineshapes of other 
rovibrational lines allow us to constrain the location of 
the corresponding 5' -wave node, and thus to determine 
the scattering length. 

The transitions which we use are from two colliding 
Na 3^S(/'= 1) atoms to the Na2 07 "purely long range" 
molecular state which asymptotically correlates to a 3^S 
and a 3^P3/2 atom [11-15]. The wavefunctions of the 
lowest vibrational levels in this potential are localized at 
distances between 50 Oo and 100 Oo, as shown in Fig. 1. 
(The Bohr radius ao = 0.0529177 nm.) This molecular 
potential is determined almost entirely by the known 



Na+ 




R(ao) 



Fig. 1. Sketch of the two-step photoassociation/molecular ionization 
process used to obtain the spectram of the 07 state. Two colliding 
atoms approach along the ground state molecular potential and are 
excited to a bound molecular state by a laser photon (solid arrow). The 
excited molecules thus created are then excited to an autoionizing 
continuum by the second laser (dashed arrow). The p-wave ground 
state scattering wavefunction is shown with a node directly underneath 
the Re of the Og potential. This leads to an absence of p -wave features 
(odd rotational lines) in the experimental spectrum of the v = vibra- 
tional level. 



long range forces between atoms and the magnitude of 
the atomic spin-orbit splitting, and thus may be calcu- 
lated to high precision. The transition rate depends on 
the overlap between the ground state wavefunction for a 
low energy collision and the excited bound state wave- 
function. It is a fortuitous coincidence that there is a 
node in the /j-wave scattering wavefunction that is 
nearly centered on the minimum of the OJ potential. 
This leads to an almost complete cancellation of the 
overlap integral between the Na ^S(/'= 1, m= — 1) + 
Na ^S(/= 1, OT= — 1) p-wave scattering wavefunction 
and the symmetric v = vibrational wavefunction, re- 
sulting in a striking and characteristic absence of 
p -wave features in the spectrum of the v = level of the 
Og" state in our experiments. We are able to construct a 
family of ground state potentials consistent with the 
known spectroscopy of the molecular ground states that 
also reproduce the /j-wave node near the minimum of 
the 07 state. We obtain further constraints on the accept- 
able potentials from the width and the relative heights of 
the rotational features in the spectrum. This, in turn, 
places constraints on the position of the corresponding 
5' -wave node. Finally, we relate the 5' -wave nodal posi- 
tion to the scattering length. 



2. Experimental Spectra and Lineshapes 

The experiments are performed by loading Na atoms 
into a "dark spot" MOT [16]. The trapping lasers are 
turned off for brief periods (~ 10 (jus) and a tunable 
probe laser is introduced during this time. For selected 
frequencies of the probe laser, red of the atomic reso- 
nance, pairs of atoms undergoing collisions are excited 
to molecular states. These molecules are then detected 
by ionization with a second probe laser. The ionization 
laser is tuned to be non-resonant with any photoassoci- 
ating transition but to allow ionization of the molecular 
states of interest. Measurements such as these have been 
described before [8,15], and here we review only those 
features important for the understanding of the analysis 
below. 

The MOT captures Na atoms using the 
3^S(f= 2) -> 3^Fy2(f= 3) atomic transition. This transi- 
tion is not a closed cycling transition because occasion- 
ally atoms get excited to the 3^P3/2(/= 2) state which can 
decay to the 3^S(f= 1) state, requiring the "repumping" 
of atoms that fall into the 3^S(f= 1) ground state. The 
dark spot MOT has this repumping frequency missing 
from the central volume of the trap and, consequently, 
the atoms are almost completely optically pumped into 
the 3^S(f= 1) ground state. All of the transitions we 
discuss in this paper begin from the 3^S(f= 1) + 
3^S(f= 1) ground state. When the photoassociating 
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probe is introduced there are no excited state atoms 
present. The ionizing laser present during the probe 
periods is tuned blue of the atomic resonance frequency 
and does not affect the atoms in the MOT. The ionizing 
laser frequency is chosen and kept fixed while the pho- 
toassociating laser is scanned over the ~ 1 GHz fre- 
quency range spanned by the rotational structure of a 
given Og" vibrational level. We check that the laser pow- 
ers are low enough that the signal heights are linear and 
that the linewidths are independent of power. 

The frequency of the ionizing laser is chosen to take 
the molecules formed in the photoassociation step into 
the ionization continuum (see Fig. 1) just above the 
3^P3/2 + 3^P3/2 asymptote. This continuum has structure 
[8] which complicates the interpretation of the spectra 
presented here. If the sum of the two laser frequencies 
(photoassociating plus ionizing) coincides with a nar- 
row feature in the continuum for some particular fre- 
quency range of the photoassociating laser then the rel- 
ative intensities of the rotational lines will not be 
proportional to the transition strengths in the photoasso- 
ciation step. Since these relative transition strengths are 
important for our analysis, we work in a region where 
there are no sharp resonances and the ionization contin- 
uum is not rapidly varying. Nonetheless, this does lead 



to some uncertainty in the relative intensities of the 
experimental peaks. 

Figure 2 shows spectra of several OJ vibrational lev- 
els. Several observations can immediately be made. The 
spectra show a rotational progression of lines at posi- 
tions given by B^J'iJ' + 1), where only the lowest five /' 
features are visible (/' = — 4), and B^ is the rotational 
constant for vibrational level v . The J' = 2 peak is always 
much larger than the other rotational lines. For the v = 
vibrational level the odd J's are nearly absent, while for 

V = 1 these odd /' peaks are clearly visible. In fact the 
odd /' peaks are larger than the /' = and 4 lines. The 

V = 5 spectrum is typical for the v > 2 levels. Moreover, 
for V = the ratio of the heights of the /' = 4 and the 
J' = 2 peaks is of the order of 0.2. Changing the fre- 
quency of the ionizing laser can change this ratio by 
approximately a factor of two. Finally, for all the vibra- 
tional levels examined up to v = 8 the J' = 2 peak, with 
a width of ~ 30 MHz, is narrower than the /' = 4 peak 
and is more symmetric as well. 

The observed lineshapes are understood as a 
Lorentzian profile convolved with the thermal distribu- 
tion of the ground state collision energies [6]. The 
lineshape for a given vibrational-rotational level (v, /') 
is proportional to the following lineshape factor: 




-200 200 400 -50 50 

Photoassociation Laser Frequency [MHz] Detuning from ttireshold [MHz] 

Fig. 2. Experimental rotational progressions for the v = 0, 1 , and 5 vibrational levels 
of the Og state. Each panel spans 900 MHz except for the lower right which is an 
expanded comparison of the v = and v = 1 , 7 = 2 peaks, showing that the v = 
feature has a larger width. The fitted curves are i-wave [see Eq. (2)] for 7 = and 2, 
p-wave for 7 = 1 and 3, and J-wave for 7 = 4, except for the v = 0, 7 = 2 peak for 
which there is a strong J-wave contribution. The temperature is fixed at k^Tlh = 9 
MHz and the natural line width is set to 20 MHz for v = and 22 MHz for v = 1 and 
5 (to allow for unresolved hyperfine structure). 
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S{oj,T,v,J')= 2 ««(2F'+1) 

F'p'^,Fptfa 






where co is the laser frequency, T is the temperature of 
the sample, Efyp, Icfifyp), and jy are the excited state 
energy, wavefunction, and natural linewidth respectively. 
The excited state wavefunction is labeled by the total 
angular momentum quantum number F', the parity p' 
and the remaining hyperfine and electronic degrees of 
freedom labeled jS. In addition, it is labeled with the 
vibrational quantum number v and rotational quantum 
number J', where J' = F' — I and / is the total nuclear 
spin angular momentum quantum number [13]. The 
summation over F'p'jS in Eq. (1) for a (v,J') level is due 
to the (unresolved) hyperfine structure of the OJ state. 
The ground collisional wavefunction represented by 
I ^Ip^a ) is energy normalized, the subscripts denote the 
spin channel \Fpifa) in which the collision starts, and 
the + indicates the proper scattering boundary condi- 
tions [17]. F is the ground state total angular momen- 
tum, p is the parity, and E is the asymptotic kinetic 
energy. The total angular momentum of the system can 
be written as F = /+/a +/b = ^ +f, where/, and/b are 
the asymptotic total angular momenta of the two atoms, 
^is the mechanical rotation,/ — the vector sum of /a and 
/b — is a generalized spin label, and a uniquely labels the 
remaining degrees of freedom of the asymptotic atomic 
scattering states for the 3^S(/,= l) + 3^S(/b= 1) colli- 
sion. The quantity ria is the population of the collision 
channel labeled by a. To avoid confusion between the 
atomic and molecular labels we will hereafter label indi- 
vidual atomic hyperfine states by /a or/b while /will be 
used solely to denote the vector sum of/ and/. Finally, 
i^Fp^a is the electronic optical transition matrix element 
between the ground state labeled by Fpffa and the ex- 
cited state labeled by F'p'jS. The rate yj^ is the rate at 
which the excited vJ' level produces observable prod- 
ucts, in this case, the photoionization rate by the second 
laser. Here the photoionization contributes negligibly to 
the total width: y,, <K jy. 

We assume that the absorption of the second photon 
does not modify the shape of the spectra. From chang- 
ing the color of the second photon we have seen that this 
is not always a valid assumption. Nevertheless, the mea- 
surements indicate that, over a large range of frequen- 
cies of the second laser, the relative intensities of the 
main features that we are concerned with in the spectra 
are insensitive to this. 

For ultracold atom-atom collisions the matrix element 
of the dipole moment has a kinetic energy dependence 
governed by the Wigner-threshold law [18,19], that is. 



the initial collision wavefunction I '^Fp'efa) is proportional 
to £'<^''*^'''*. For example, for ^ -wave scattering the wave- 
function is proportional to vF. Due to this Wigner-law 
variation in the (Franck-Condon) matrix element, Eq. 
(1) leads to asymmetric lineshapes [6] where the blue 
side is dominated by the Lorentzian in Eq. (1) and the 
red side is predominantly determined by the Maxwell- 
Boltzmann distribution of kinetic energies. The ob- 
served position of the peak is always red shifted with 
respect to the actual bound state energy Eyr. This shift 
is on the order of hgT, the linewidth is on the order of 
^bF+ Jy, and both increase with f. 
For each J' we fit the line to 



Snt{co,T,v,J')=Ayr 

r" p(2fti)/2 

ciFe"'^*B7- £ 

J„ {E + ^co - Ey,rf + {yy/2f ■ 



(2) 



The coefficient Ayr is the overall amplitude, F^y is the 
transition threshold energy, y.. is the linewidth and T the 
temperature. The results of our fits are shown in Fig. 2. 
We use a single value of F for all of the data, determined 
from the fits to be (450 ± 50) |juK (kBT/h = 9 MHz). For 
reasons discussed below, we fit the odd J' features to Eq. 
(2) with f= I (p -wave) only. The J' = and 2 peaks are 
fit to /= (5 -wave), except for v = where we find it 
necessary to use a sum of f= and f= 2 contributions. 
The J' = 4 peak is fitted with just f= 2 (d-wave). The 
natural linewidth of the OJ states is 20 MHz, which is 
twice the atomic linewidth [20,21]. For v = we expect 
the unresolved hyperfine structure to broaden the line by 
= 2 MHz. To fit the v = 0, J' = 2 peak with a single 
5 -wave lineshape requires an unrealistically large (30 
MHz) linewidth, whereas for v = 1 , where the hyperfine 
splitting is slightly larger, a linewidth of only 22 MHz is 
required to fit the data. We return to these points in Sees. 
3 and 4. 



3. General Theory 

The theory which underlies our calculation of the 
spectrum involves three major pieces: the ground state 
wavefunctions, the excited state wavefunctions and the 
molecular Rabi matrix which gives the optical coupling 
between them. These determine the transition amplitude 
matrix element {(i)Fy^\fifiFp^a\''i^Fp*?fa), from which we 
calculate synthetic spectra to compare to experiment. 

The first piece is the ground state wavefunction 
I ^i>^a ), which is obtained from an exact solution of the 
Schrodinger equation for the ground state Hamiltonian 
Aground for a given set of adiabatic Born-Oppenheimer 
(ABO) potentials which are derived from experimental 
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Rydberg-Klein-Rees (RKR) potentials. The ground state 
Hamiltonian //gfound is set up for a given value of the total 
angular momentum and parity and includes electrostatic 
interactions V{R) (the adiabatic Born-Oppenheimer po- 
tentials), the mechanical rotation operator ^llfjbR^, the 
radial kinetic energy operator, the spin-spin dipole inter- 
action, and the atomic hyperfine Hamiltonians. Most of 
our discussion will use a simpler model of //gmund and 
I '^FpfD since this provides greatly improved insight. We 
note that although the discussions may be based upon 
simpler, intuitive models the final calculations use the 
Ml H',S.uni and mi,Vf a). 

The next piece of the theory required to model the 
photoassociation spectra is to calculate the excited rovi- 
brational-hyperfine wavefunctions \4>Fyfi) and energies 
Epy^. Once again, these are obtained from an exact 
treatment of the excited state Hamiltonian HSi^d which 
includes the same interactions for the excited state as 
were contained in //pound plus a spin-orbit interaction 
that results from the presence of the excited Na 3^P 
atom, and retardation of the excited resonance dipole 
interaction. A discussion of //eSed and methods for find- 
ing its bound state solutions are found in Refs. [13] and 
[14]. Once again, most of our discussion will be based 
on a simple one channel adiabatic picture of the OJ 
bound states although the exact bound state wavefunc- 
tions and energies are used in the calculations. 

Finally, we need the molecular Rabi matrix elements 
■i^f)?^ between the initial ground electronic state labeled 
by ffa and the excited electronic state labeled by 13. 
Dipole selection rules require that p'=—p, and 
AF = F'- F=[0,±l}, except that AF ^ for F = 0. 
The Opp^^a are calculated from the known atomic transi- 
tion dipole moment between a ground Na 3^S atom and 
an excited 3^P atom using the basic approach described 
in Ref [21] but generalized here to include hyperfine 
structure. The molecular Rabi matrix elements depend 
on the excited rovibrational-hyperfine state quantum 
numbers, F'p'fSvJ', and the ground state hyperfine levels 
/a and/b of the two colliding atoms. 

These three pieces of theory are integrated together 
using Eq. (1) to yield a theoretical spectra which can be 
compared to the experimental spectra. We know that we 
can calculate the excited state 07 bound state energies to 
an accuracy of a few MHz [13] and have used this 
capability to determine a precision value of the Na 3^P3/2 
lifetime and to provide the first experimental verifica- 
tion of retardation of the interaction between two atoms 
[14]. 

Below we will briefly describe each of these three 
theoretical parts while emphasizing those portions rele- 
vant to the current problem of extracting ground state 
scattering lengths. Many arguments will take advantage 
of simple physical pictures. These pictures are meant to 



be intuitive and they have been verified within the con- 
text of two colliding Na atoms where possible. However, 
we note that the final results are based on the full 
Hamiltonian, the exactly calculated ground and excited 
state wavefunctions, and the hyperfine labeled electronic 
transition dipole moment between the initial and final 
hyperfine labeled electronic states. 

3.1 Ground State Dynamics 

Although we have set up a complete quantum scatter- 
ing calculation for two ground state atoms with hyper- 
fine structure, as described in the previous section, a 
sufficiently accurate model of ^S + ^S collisions is ob- 
tained with the atomic hyperfine Hamiltonian for each 
atom, the ground Z'2g and a^2u molecular potentials, 
the mechanical rotational kinetic energy, and the 
—ft^llpu • d^/dR^ radial kinetic energy (where the re- 
duced mass jx equals half the atomic ^^Na mass). This 
approximate model ignores the very weak magnetic 
spin-spin interactions and the second-order spin-orbit 
interaction with distant electronic states. In the absence 
of these weak spin-dependent terms in the Hamiltonian, 
the mechanical rotation ^ is a conserved quantum num- 
ber. This does not imply that ^-changing collisions are 
always irrelevant. In fact, in experiments aiming at Bose 
condensation, atom loss is in a large part due to such 
processes, which can always be treated using a weak 
interaction picture [22,23]. However, spin interactions 
play a negligible role in the description of the spectra 
obtained with photoassociative spectroscopy. 

The electrostatic Z'2g and a^2u potentials over part 
of the range of their attractive wells have been derived 
from conventional spectroscopy [24]. We extrapolate 
these RKR potentials by joining them smoothly to the 
familiar long-range dispersion form Vjisp = — 2"=6 CJR" 
using the coefficients of Ref. [25]. Note that for R> 
30 flo these two adiabatic Born-Oppenheimer potentials 
are essentially identical and are, at 30 Oo, about 
Vdisp/^B = ~ 0.7 K deep. These potentials predict that 
the Z'2g state has 65 ^-wave vibrational levels while the 
a^2u potential has 15 5- wave levels [24,26]. The scatter- 
ing length associated with each potential is sensitive to 
the precise phase of the wavefunction at zero energy, 
which is related to the binding energy of the last bound 
state. Uncertainty in the extrapolation of the RKR re- 
gion of the potential leads to uncertainty in the exact 
position of the last ground state vibrational level, and 
consequently uncertainty in the scattering length. It is 
the sensitivity of the photoassociation spectra to the 
phase of the low energy ground state wavefunction (i.e., 
to the position of the nodes in the wavefunction) that 
allows us to obtain the scattering lengths associated with 
the collision of particular hyperfine states. In order to 
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reproduce the experimental 0~ lineshapes we will allow 
the shape of the inner wall of the electrostatic Z'2g and 
a^2u potentials to vary in order to adjust for short and 
long range extrapolation uncertainties, but we restrict 
the changes to conserve the number of levels in these 
two potentials. In practice, the inner walls of the two 
RKR curves are allowed to vary independently. 

In the dark spot MOT the sodium atoms are in the 
atomic /a = 1 hyperfine state and are assumed to be 
distributed equally over the three magnetic sublevels m,/^. 
Since the MOT has a nearly-zero magnetic field (< 0. 1 
mT and spatially-varying in magnitude and direction), 
collisions are independent of the orientation of the 
molecule in the laboratory frame. We may view the 
collision as starting when the atoms are infinitely far 
apart with a definite value for the relative angular mo- 
mentum f and retaining this value throughout the colli- 
sion. We can therefore evaluate the ground state Hamil- 
tonian in the atomic hyperfine basis \Fpifa) for fixed 
values of the total angular momentum F = i+f and 
parity p , where here a designates {//b } . The parity p is 
the symmetry of the ^S + ^S Hamiltonian under inver- 
sion through the center of mass of all the electron and 
nuclear coordinates. Since the angular momentum i is 
conserved during the collision, coupling to F is not 
really necessary but is useful in setting up the molecular 
Rabi matrix below. The rotational and hyperfine Hamil- 
tonian terms are diagonal in this atomic hyperfine basis, 
although the electrostatic terms are not, since the basis 
does not form states with good electron spin S =s^ + s^- 
However, when we neglect the weak spin- spin coupling 
terms, there is a diagonal representation in a molecular 
basis with S and i as good quantum numbers: 



\Fp£fSI) oc 2 V(25+l)(2/+l)(2/,+ l)(2/i,+ l) 



//b 



X. 



''a 'a J a. I 

% ib /b \\' 



Sb lb fb \\FpifMb) 
S I f] 



(3) 



where {...} is a nine-j symbol; the exact equation has 
phase and normalization factors resulting from nuclear 
symmetrization. Since the Born-Oppenheimer curves 
do not depend on / it is a conserved quantity. There is 
also a restriction on the permissible quantum numbers 
due to the homonuclear nature of the dimer since the 
basis states must be antisymmetric with respect to ex- 
change of the two nuclei. This leads to the restriction 
(— 1)''+'^+' = 1 with o" = 0(1) for gerade (ungerade) states 
(for ^S + ^S collisions there also exists a one-to-one cor- 
respondence between gerade/ungerade and the total 



electron spin S, allowing S to be substituted for cr). In 
the atomic basis the restriction is (— l)'*^"^-"^'' = 1. 
An important consequence is that the Na^S(/a = 1) + 
Na^S(/b =1) spin state couples to even/= or 2 for even 
partial waves and to odd /= 1 for odd /"s. This latter 
statement is true whether or not we neglect the weak 
spin-spin interactions. 

The fact that f and / are good approximate quantum 
numbers lets us develop a relatively simple picture of 
photoassociation spectra due to collisions of ^S(/a = 1) + 
^S(/b =1) atoms. There are only two possible 5 -wave 
contributions, corresponding to /= and /= 2. These 
have F = and F = 2 respectively. For the p -wave there 
is only one possible contribution, corresponding to/= 1 
and F=0, 1, or 2. Finally, there are two possible cf-wave 
contributions, where F = 2 for/= and F = 0, 1 , 2, 3, or 
4 for/= 2. Within our approximation of neglecting weak 
spin-spin interactions a given/, /subspace contained in 
Hamiltonians labeled by different F's are identical, with 
identical wavefunctions. Thus, the three values of F 
which contain the/= \, i=\ subspace have identical 
;? -waves and thus identical nodes. Therefore, we can 
represent the collision in terms of two 5 -waves, one 
p-wave, and two cf -waves. For brevity we will refer to 
these five wavefunctions as ^|/' and thus as ^J\ ^5', 
^pV, ^.t', and ^^>. 

BEC experiments can magnetically trap the alkali- 
metal atoms in one of the magnetic sublevels. There are 
two relevant states. One is the doubly polarized state 
where all atoms are in the atomic /a = 2 and nt/^ = 2 state. 
Two of these atoms have a projection of mf=mf^ + 
m/^ = 4 which implies /= 4. The second trappable spin 
state, used by the MIT group [2], is the /a= 1 and 
ntf^ = — 1 state. This implies that a collision between 
two such states couple to a I (/a = 1 /b = i)f='2;mf= — 2) 
state. The zero-field scattering length of the latter state 
is extracted from our experiment; in fact, it is related to 
^5'. Because the magnetic fields used in the sodium 
traps of Ref. [2] are weak, the Zeeman shifts of the 
atomic hyperfine states are small compared to the hy- 
perfine structure and thus have little effect on the colli- 
sion dynamics. Hence the zero-field scattering length is 
the relevant parameter in those experiments. 

The ^S + ^S collisional wavefunction is inherently 
multichannel. In Fig. 3 we show the three components 
of an exact close-coupling wavefunction [22,23,27], for 
an incoming s -wave in the /a = 1 , /b = 1 channel with 
/= F = 2 and a kinetic energy of E/ks = 500 (xK. The 
figure also shows the three potential curves (dashed 
lines) for each of the three spin channels. The horizontal 
line indicates the total collision energy. The plane wave 
scatters into the two other 5 -waves with/= F = 2; they 
have /a = 1 , /b = 2 and /a = 2, /b = 2 respectively. These 
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Fig. 3. The multichannel S + S collisional wavefunction as a function 
of internuclear separation. The wavefunction describes an E/ks = 500 
(jlK i-wave collision of two /« = 1 atoms coupled to a /=/« +/b = 2 
state. Two of the three spin components ip of the wavefunction decay 
exponentially because those states are asymptotically unaccessible. 
The dashed lines denote the attractive long-range dispersion potential 
for each of the three spin channels. The horizontal line denotes the 
total energy in the collision. 

Other channels are closed asymptotically by E/ks = + 85 
mK and +170 mK, respectively. Therefore, they are 
only populated at short internuclear separation, where 
the attractive potential is larger than the asymptotic sep- 
aration and where the electrostatic exchange interaction 
(the difference between the Z'Sg and a^Su potentials) 
can mix these three spin channels. The mixing occurs 
around 25 Aq, where the exchange splitting is compara- 
ble to the hyperfine splitting. Inside 20 Oo the wavefunc- 
tion oscillates rapidly due to the high kinetic energy in 
the deep potentials and shows striking interference pat- 
terns due to the strong electrostatic interaction. In this 
region the ' 'molecular' ' basis would be more appropri- 
ate than the atomic hyperfine one. For R>30 Oo the 
three channels are decoupled and the dynamics is gov- 
erned by the common long-range potential and the ki- 
netic energy. The wavefunction components for the up- 
per two channels decay to zero since these channels are 
closed, while the 5' -wave in the /a = 1 +fi, = 1 entrance 
channel extends to R = °° with long wavelength oscilla- 
tions. At large R this low-energy wavefunction (except 
for an R independent phase factor) is given by 



-^ 



Ifji 1 



^"ttV^ 



sm(k(R — fli.-i)) 



\F=2,p = +\,-e=0J=2(f,= \,f^=\)), 7?^»(4) 

with k the asymptotic wavenumber and Oi-i the scatter- 
ing length. 

Most notable about the wavefunction in Fig. 3 is the 
node around 60 Oo and the absence of appreciable prob- 
ability in the two asymptotically closed channels for 
internuclear separations larger than 50 Oo- In the rest of 
this paper we adopt the convention of calling this node 
the last node in the wavefunction, even though the wave- 
function keeps oscillating with a wavelength corre- 
sponding to a kinetic energy of 500 (xK. The E = 
wavefunction will always have a last node associated 
with the number of bound states in the potential (see 
Appendix A), and this nodal position does not change 
significantly for wavefunctions with kinetic energies be- 
low 1 mK. A more general expression for the asymptotic 
wavefunction in Eq. (4) replaces sin(k(R — ai,-i)) with 
sin(kR + 8(k)) where the phase shift 8 has as a limiting 
behaviour — ai-,k for small collision energies. The an- 
swer to the question ' 'what is small?' ' is system-depen- 
dent, but for Na the answer is about 1 mK or less. 
Moreover, for these collision energies and for internu- 
clear separations R at which the long-range dispersion 
potential has died off sufficiently compared to the ki- 
netic energy, the product kR is still small compared to 
one and the wavefunction in Eq. (4) can be approxi- 
mated as being proportional to k"^(R — Oi- 1). The 
wavefunction for higher-order plane waves is propor- 
tional to A:'^''*^''^. This analytic variation with k defines 
the Wigner threshold regime [18,19]. 

In Fig. 4 we show the radial density of three ground 
state wavefunctions as a function of internuclear separa- 
tion. All wavefunctions correspond with a collision start- 
ing in a/a = 1, /b = 1 channel with 500 (xK kinetic en- 
ergy. The density is obtained from the multichannel 
wavefunction I '^Fp'^k^^) ^Y summing the squares of the 
4'fJb (^ ) ^t each R . In particular, the graph shows the 
^2', ^1\ and '^di waves. Moreover, Fig. 4 shows the 5'-, 
p-, and flf-wave potentials of the /^ = 1 ,/b = 1 component 
of the potential matrix. In the radial region that is impor- 
tant for the photoassociation spectroscopy of the 07 
state this diagonal element of the multichannel potential 
matrix is given by - Ce/R'' + (;?V2^)^(^+ l)/R\ This 
is a consequence of the fact that for these internuclear 
separations the two ABO potentials are identical and 
given by their dispersion form. Moreover, the density for 
R> 50 ao is solely due to/a = 1,/b = 1 component of the 
wavefunction. 
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Fig. 4. a) The s,p, and J-wave potential barriers as a function of internuclear separa- 
tion, b) The probability densities for the three wavefunctions, corresponding to a 500 
H,K collision starting from the /„ = 1, /b = 1 spin channel. 



For Na the height of the fi?-wave barrier maximum at 
75 flo is 5.4 mK. This is much higher than the tempera- 
ture (~ 500 |jtK) of the atoms in the MOT. Therefore, 
the penetration of the fi?-wave into the region near 75 Oo 
is greatly reduced by the centrifugal barrier. In fact, full 
close-coupled calculations show that, for Na MOT tem- 
peratures, £> 1-wave wavefunctions outside of the bar- 
rier are almost independent of the shape of the electro- 
static potentials inside the centrifugal barrier. Therefore, 
the fi?-wave wavefunction is mainly determined by the 
well-known long-range form of the potential while 
higher partial waves do not contribute significantly to 
the lineshapes. As a result, we find that "^al and also 
^o' are almost identical to a purej^C^^) spherical Bes- 
sel function in the region where the Franck-Condon 
factors are nonzero (i.e., in the region of the centrifugal 
barrier) with their normalization determined by asymp- 
totic boundary conditions. This implies that we will have 
no freedom in modifying the fi?-wave features of the 
spectra. There is much more penetration of the s- and 
p-v^ssie. wavefunctions to small internuclear separations 
and therefore they will display a significant dependence 
on the shape of the inner wall of the two ABO poten- 
tials. 



The above is in contrast to the case of ^'Rb where a 
fi?-wave shape resonance dominates the spectrum ob- 
tained from samples of doubly polarized atoms [28]. In 
Rb, the fi?-wave barrier is comparable to the most proba- 
ble collision energy (k^T) and as a result there is signif- 
icant barrier penetration by the wavefunction. A similar 
effect could occur in the current Na experiments for the 
/7-wave; however, this is in contradiction with the obser- 
vation of & p -wave node near the minimum of the OJ 
state. Because the d-wave barrier height in Na is large 
compared to the most probable collision energy, any 
fi?-wave resonance that might occur will be narrow. No 
experimental evidence exists for such a resonance. 

3.2 Excited Bound States 

The long-range Og" potential results from a spin-orbit 
avoided crossing between a ^2g and a ^Ilg potential [11- 
13]. These two non-relativistic electronic curves plus six 
additional potentials dissociate to the atomic ^S + ^P 
asymptote [11]. The notation ^^*'Ao^ reflects the underly- 
ing symmetries in the nonrelativistic electronic Hamil- 
tonian, for which the total electron spin S is conserved 
since the electrostatic interactions are independent of 
spin. The absolute value of the projection of the total 
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electronic orbital angular momentum on the body-fixed 
symmetry axis (A) is conserved due to the cylindrical 
symmetry of the electronic Hamiltonian. The labeling 
of the molecular states with a, which is either gerade (g) 
or ungerade (u), is a result of the inversion symmetry of 
all electrons through the center of mass of the molecule. 
Movre and Pichler [11] showed that if one constructs 
a Hamiltonian based on both electrostatic interactions 
and the relativistic spin-orbit interaction that results 
from the P atom, then the resulting Hamiltonian mixes 
electronic states labeled by SAXcr (where X is the body- 
fixed projection of S) with states labeled by S'A'X'cr' 
such that 12 = A + X= A' + X' is conserved and a= a'. 
In addition, for i7 = states the Hamiltonian also sepa- 
rates into two subspaces which have definite symmetry 
under reflection of the electronic wavefunction through 
an arbitrary plane containing the internuclear axis. This 
reflection symmetry is denoted by a superscript + or — . 
The complete notation for the spin-orbit mixed Hund's 
case(c) states is 12^. The purely long range Og potential 
is obtained within this two-state Movre-Pichler model 
by incorporating only the spin orbit and resonant dipole 
interactions which are the dominant forces at long range 
between an alkali ^S atom and a ^P atom. The two 
adiabatic OJ potentials are found by diagonalizing the 
potential matrix: 



Vm 




2 
Via . 

2C3 A 
R^ 3 ■ 



n 



(5) 



where A is the atomic spin-orbit splitting and we have 
taken the zero of energy to be the ^S + ^¥3/2 asymptote. 
Within this simple model the well depth is A/9, indepen- 
dent of the resonant dipole interaction strength and the 
potential minimum is at R^ = (9C3/2/I)'". For Na(3^ P), 
A = 515.520 GHz, C3 = 4.018 zJ nm' (6.219 a.u.) [14] 
and Re ~ 72 Aq. 

Figure 5 shows the purely long range adiabatic 0~ 
potential along with the three lowest adiabatic vibra- 
tional wavefunctions in this potential. This is a purely 
long range potential in the sense that the electron clouds 
of the two atoms do not overlap in the vicinity of the 
potential well and it is therefore completely determined 
by atomic parameters. In the region where these wave- 
functions are nonzero, the OJ potential is nearly a har- 
monic potential and hence, the v = and 2 wavefunc- 
tions are nearly symmetric with respect to R^ while v = 1 
is antisymmetric. 

In Ref. [13], three of the present authors discussed the 
rotational and hyperfine structure of the OJ vibrational 
levels. There, we showed that we could obtain the exact 




R (units of a ) 



Fig. 5. The purely long-range Og adiabatic potential and selected 
adiabatic vibrational wavefunctions versus the internuclear separation. 
The V = and 2 wavefunctions are nearly symmetric with respect to 
the minimum of the well, while the v = 1 level is antisymmetric. 

bound states of the fully rotating 3^S + 3^P Hamiltonian 
including hyperfine structure. For a given total angular 
momentum F' and parity;?', the full Hamiltonian matrix 
will include up to 96 coupled-spin basis states. Although 
the multichannel wavefunctions in principle can be dis- 
tributed over as many as 96 spin channels, an appropri- 
ate transformation can usually be found that will con- 
strain the nonzero amplitude to at most a few channels. 
Moreover, the nonzero components of such a wavefunc- 
tion have a common radial dependence, as depicted in 
Fig. 5. In other words, the Og levels for v < 9 are essen- 
tially adiabatic [13] and thus can be viewed as single- 
channel wavefunctions. Note, that the actual spin struc- 
ture is essential for calculating the transition matrix 
elements which are labeled by F'p'jS. 

For the purely long range Og states, it turns out that 
the hyperfine and Coriolis interactions are absent in first 
order. Therefore, in addition to the quantum numbers F' 
and p', the quantity J' = F' — I is approximately good. 
Moreover, J' = S + L+ €, where the electron orbital an- 
gular momentum L=\. General symmetry relations 
show that p'={—iY*^ for homonuclear ^S + ^P 
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molecules, i.e., odd { corresponds to even parity and 
vice versa, and (— 1)''+°^+' = — 1 where a = 0(1) for ger- 
ade (ungerade) states. For the OJ states where S = L= 1 
additional selection rules are appropriate. We find /' + / 
is odd or, equivalently, even /' correspond with odd 
parity and vice versa. Some of these selection rules 
slowly break down with increasing vibrational quantum 
number as second order coupling to nearby states with 
different O^ symmetry becomes stronger. 

This description of the Og" vibrational levels leads to 
the following picture of the level structure. The energy 
level distribution is in first order given by a rotational 
progression in /'. Each /' consists of a group of nearly 
degenerate levels. The /' = level is two-fold degenerate 
with / = 1 or 3, while the /' = 1, 2, 3, and 4 levels are 4, 
8, 6, and 10-fold degenerate, respectively. From Ref. 
[13] we know that for the lowest three vibrational levels 
the hyperfine degeneracy is lifted by no more than 5 
MHz, which is still small compared with the natural 
width and the rotational constant. 

Even though /' is a good approximate quantum num- 
ber and behaves as an effective rotation, this does not 
imply that states with a definite value of the mechanical 
rotation are formed. In fact, even J"s represent positive 
parity states and therefore contain even partial waves 
and odd /"s contain odd partial waves. For example a 
/' = 2 state will have /= 0, 2, and 4 contributions. The 
low temperatures in the present experiments limit f to 
values of 2 or less. 

3.3 Molecular Rabi Matrix 

The molecular Rabi matrix elements /2|^*'^ are ob- 
tained by first considering the allowed optical excitation 
of a pair of atoms by a single photon at large internu- 
clear separation. The Rabi matrix in the atomic hyper- 
fine basis is then transformed into the molecular basis. 
The basic approach is an extension of that originally 
used in Ref. [21 ] where we have incorporated the atomic 
hyperfine structure. In simple terms, we know the 
atomic transition dipole moment and the atomic hyper- 
fine selection rules for optical transitions, which are 
A/a = {0, ± 1 }, AZa = 1, Ai'a = 0, and A4 = 0, where we 
have assumed that the atom labeled "a" has been ex- 
cited. These selection rules insure that only the orbital 
angular momentum k changes for the optically dipole 
allowed ^S ^ ^ P transition. 

At large internuclear separation we can define a set of 
atomic scattering states 



l/^,c/am/^,c/b'n/i) = Yeu,\c/^mf)\c^ft^f^) 



(6) 



chanical rotation of the two atoms about their 
center-of-mass. In the above description c^ stands for all 
other quantum labels needed to uniquely specify the 
atomic hyperfine state — i.e., for a ground state Na atom 
Ca = 3^S while for the first excited state of Na c^ = 3^R 

Beginning with an initial set of atomic scattering 
states l//x,Ca = 3^Sfjnf^, c^ = 3^S/bm/^) and a second set 
of atomic scattering states l/'/x',Ca = i^V flm}^, c{, = 3^S 
fbrtif^), where we arbitrarily assume that atom "a" is 
excited, then it is obvious that we can derive the Rabi 
matrix elements between these two states from the 
known atomic transition dipole. In such a picture 
the Rabi matrix element will be zero unless 
S//S^b,cbS/^/i,S„^i_,„yi_= 1, and the hyperfine selection 
rules for the optically excited a -atom are obeyed. These 
selection rules insure that only one atom absorbs the 
photon when the two atoms are at infinite internuclear 
separation. The real situation is slightly more compli- 
cated since we must symmetrize the asymptotic basis 
with respect to exchange of the identical nuclei. 

Our asymptotic derivation of the molecular Rabi ma- 
trix is strictly valid for the purely long range Og state, 
since the electronic clouds of the two atoms never over- 
lap and distort the atomic dipoles. As a check on the 
transition dipole moment and a confirmation of our code 
we can calculate the natural lifetime of an arbitrary 
molecular state; e.g., the A '2u state or the purely long 
range OJ state. This involves summing over all ground 
state hyperfine components and, as expected, yields 
~ 20 MHz for the purely long range Og" state and — 10 
MHz for the A 'Xt state. 

3.4 Evaluation of the Molecular Transition 
Strength 

The absorption of a photon excites the colliding 
atoms from a ground state scattering wave into a bound 
excited state molecule. Although our analysis is based on 
exact numerical calculations of the molecular Rabi ma- 
trix and the ground and excited state multichannel quan- 
tum wavefunctions, much physical insight for interpret- 
ing our result can be obtained from considering the 
molecular transition strength labeled by the approxi- 
mately good quantum numbers discussed above: /', /, 
and /. This transition strength is determined from the 
Franck-Condon overlap matrix elements: 



F}f{E)= 2 \{4>Fvii\fiOf;imX)\^ 



(7) 



a^F'F 



which are products of magnetically resolved atomic hy- 
perfine states kj'jnfj) for atoms a= {a,b} and a spher- 
ical harmonic wavefunction F/^ which describes the me- 



The sum over a only involves channels where the two 
atoms have /a = /b = 1 . The summations over p and p ' are 
absent as f uniquely defines the parity of the ground 
state and p' = — p from the selection rules of the transi- 
tion dipole moment. 
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The discussion in Sec. 3.1, when combined with the 
above equation, shows that there are only two possible 
5 -wave contributions, corresponding tof=F=0 and 
f=F=2. These are designated as F^o(E) and Fli(E), 
respectively. For the purely long range OJ state, the 
5 -waves contributes predominantly to the J' =2 and to a 
lesser extent to the /' = feature. For the /j-wave there 
is only one possible contribution, F'pi{E). The p-wave 
contributes to /' = 1 and 3 features only. Finally, there 
are two possible d-wave contributions, f/oiE) and 
Fd2(E). The fi? -waves contribute to the /' = 0, 2, and 4 
features. 

One important aspect of our argument below is that 
fAXE) » fJiiE). Therefore, the analysis of the 
lineshapes is primarily sensitive to the/= 2 5-wave and 
not the/= one. One reason for this is that the phase 
space factor 2F' + 1 is much larger for the/= 2 5 -wave. 
However, there is no reason why the scattering length 
O/^o should be the same as the scattering length a^2, 
since the different / values lead to slightly different 
Hamiltonians. Both of these scattering lengths are dif- 
ferent from those for the electrostatic potentials for the 
'2g and ^2u states without hyperfine structure, because 
of the strong mixing of these states in the 5 -wave colli- 
sion for a given/. Our complete close coupling calcula- 
tions show: 1) that a^ is actually near 0^2, crossing it as 
the inner ABO potentials are varied, and 2) that 
Fl2(E) » Fli(E) is valid for the transitions we study. 

Finally, we make a more quantitative argument that 
near R^ the harmonic nature of the OJ potential for 
V = — 2 (Fig. 5) helps explain the relative intensities of 
the;? -wave features for these levels. Consider the follow- 
ing one-dimensional spinless Franck-Condon factor: 



f 



dRcl,,(R)^;XR) 



(8) 



In this equation, cfvi^) is the adiabatic 0~ vibrational 
wavefunction and, as discussed above, ^V is the single 
p-wave for/a = 1 +/b = 1 collisions. We neglect any R- 
variation in the Rabi matrix elements for different hy- 
perfine components of the upper level. The v = func- 
tion, and to a lesser extent the v = 2 function, is nearly 
symmetric about the minimum near R^ = 72 Oo, whereas 
the v= 1 function is antisymmetric. Since the p-wave 
has a node so close to R^, it also is nearly antisymmetric 
about R^. Therefore, the molecular transition strength for 
p -waves is very small for v = and 2, but much larger 
for V = 1 . 



4. Obtaining the Scattering Length 

Having developed these theoretical tools, we now re- 
turn to the interpretation of the experimental spectra in 



terms of the s, p, and d wavefunctions. As explained in 
Sec. 3, there are three theoretical elements which are 
needed in order to simulate the experimental spectrum 
using Eq. (1). These are the ground state wavefunctions 
I^Fp^a), the excited state wavefunctions I<^f>p), and the 
molecular Rabi matrix elements /2|^*'^ . Because of the 
checks on the transition dipole moment described in 
Sec. 3.3 we can be confident in the determination of the 
latter. Refs. [13,14] on the rovibrational-hyperfine states 
of the Na2 OJ state provide compelling evidence that we 
can calculate the excited states accurately. Thus, the 
uncertainty in our ability to simulate the experimental 
spectra is mainly associated with inaccuracies of the 
Z'2g and a^Xu RKR potentials, and thus in generating 
the ground state wavefunctions. 

In Fig. 6, we show the v = simulated spectrum for 
our original fit of the ground state Na2 RKR potentials 
[24]. The ground state collision wavefunctions are com- 
puted exactly given these Z'2g and a^2u potentials. The 
three elements of the theory are then substituted into Eq. 
(1) and the thermal lineshape is calculated assuming a 
temperature T = 450 (xK. Note that unlike the experi- 
mental spectrum (Fig. 2) the simulated spectrum has 
very large / = 1 and 3 peaks and a rather weak J =2 
feature. The reason for this is that our fit of the Na2 X'2g 
and a^Xu RKR potentials caused ^2' to have a a,--, 
scattering length of 73 Ug, with a corresponding i-wave 
node at 78 Oo- This results in a nearly zero Franck-Con- 
don factor for the 5 -wave J =2 feature. For these poten- 
tials the p -wave node for ^V was at 95 Oo, far from R^. 
This is inconsistent with the experiment and indicates 
that the RKR potentials must be altered. 




-54.2 -54 

Energy(GHz) 



-53.6 -53.6 



Fig. 6. Simulated v = spectra for the original RKR potentials. The 
/= 2 peak is completely dominated by the J-wave scattering, as can 
be inferred from its relatively large width of «= 60 MHz and the 
"slow" onset of the red side of the line. 
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Changing the inner wall of the Z'2g and a^2u RKR 
potentials changes the accumulated phase of the wave- 
function or, equivalently, changes the position of the last 
node. In Fig. 7, we show how varying the inner walls of 
the potentials modifies various properties which depend 
on the ground state scattering wavef unction. The two 
axes represent independent, adjustable parameters 
which cause a smooth change in the inner wall of the 
Z'2g and a^Xt potentials respectively. The precise form 
of the adjustable parameter is irrelevant [29] since we 
are only sensitive to the accumulated phase up to the 
Franck-Condon region (R > 50 Oo), where the potentials 
are completely determined by atomic properties. The 
plotted lines forming two distinct bands correspond to 
lines of constant position of the last p -wave node and 
constant ratio of the J =2 and /= 4 peak heights. The 
intersection of the bands in Fig. 7 determines the al- 
lowed range of the scattering length. 



larger a^ ^j 




E parameter 



Fig. 7. Parameter space plot for variation of X'Sg and a^Xi potentials. 
Note that smaller values of the parameters imply a deeper potential. 
The arrow indicates the direction in which the ai -i scattering length 
increases. 

Fig. 8a shows how the simulated spectrum changes 
when the p -wave node moves to smaller R for nearly 
constant Oi- 1 scattering length. The spectra have been 
normalized with respect to the J' = 2 peak. Notice that 
a relatively small change in the p -wave node position 
has a marked effect on the odd /' peaks in the spectra. 
Hence to have very weak v = 0, /' = 1 and /' = 3 peaks, 
consistent with the experimental data, we find that ^t' 
must have a node close to ^e- The calculations strongly 
constrain the p -wave node to 73 Aq ± 3 Aq- This defines 
the ;? -wave band in Fig. 7. Note that there is a range of 
Z'2g and a^2u potentials which satisfy this constraint. 
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Fig. 8. Simulated v = 0, Og^ spectrum for various potentials. The exact 
transition dipole moment is used: a) shows the effects of moving the 
p-wave node while ai-i is held nearly constant and b) shows the 
effects of moving ui -i while the p-wave node is fixed at 73 oo- 

In the discussion of the optimal position of the last 
p-wave node we used the wavefunctions with 500 (xK 
kinetic energy in the incoming spin channel. Unlike for 
5 -wave scattering, where in the Wigner threshold regime 
the nodal positions are independent of the collision en- 
ergy, the position of the p -wave node always shifts with 
collision energy. In fact, the zero energy wavefunction 
has a node which is about 2 Oo to 3 Oo inside the reported 
p-wave node. The 500 |jiK collision energy is close to 
the most probable collision energy in a MOT, and there- 
fore the spectra are most sensitive to the position of this 
node. 

Having determined the position of the last p-wave 
node, we now argue that the corresponding /= 2 i-wave 
node lies at smaller R. This has been confirmed by 
independent full close coupled calculations, by the theo- 
retical arguments presented in appendix A, as well as 
being supported the widths of the observed lines. Ap- 
pendix A also gives an analytical one-to-one correspon- 
dence between the 5 -wave nodes and the scattering 
length. For now it is sufficient to keep in mind that for 
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Na the value of the scattering length is always a few Co 
smaller than the position of the last node. 

In Fig. 8b the simulated spectra for several trial 
ground state potentials are shown, keeping the ^t' 
p-wave node fixed. Once again, the spectra have been 
normalized with respect to the J' = 2 peak. The figure 
shows that the /' = 4 to /' = 2 peak ratio varies dramati- 
cally with the fli,- 1 scattering length. If this were the sole 
difference we could not be as confident about our final 
values since experimentally we have seen as much as a 
factor of two change in the /' = 4 to /' = 2 peak ratios by 
varying the frequency of the ionization laser. In the 
simulations, changing the ai,_i scattering length while 
keeping the /j-wave node fixed also causes a large 
change in the width of the v = /' = 2 feature. This is 
because the width is determined from a mixture of s- 
and flf-wave contributions: an increased of -wave contri- 
bution implies a larger width. The /' = 2 width de- 
creases with decreasing scattering length because the 
flf-wave contribution becomes less and less important as 
the i'-wave Franck-Condon factor increases. Thus, the 
width of the / = 2, v = feature can also be used in 
constraining the scattering length. 

As explained in Sec. 3.1, the of -wave wavefunction is 
given by a spherical Bessel function, J2(kR)/yk -^ 
k^'^R ^ as A: ^ 0, independent of the shape of the poten- 
tial because the centrifugal barrier inhibits penetration 
of the wavefunction into the region of interest, as seen in 
Fig. 4. Thus, the intensities of the of -wave features in our 
simulated spectrum are fixed. This has been confirmed 
computationally for all the various potentials used in this 
modeling. However, changing the i'-wave node and 
thereby the Ci . i scattering length changes the ampli- 
tude of the i'-wave scattering wavefunction in the vicin- 
ity of the minimum of the Og potential, and thus the 
strength of the i'- wave features. Moreover, as the i'-wave 
character of the v = 0, /' = 2 peak increases, the 
linewidth of the feature becomes narrower. Therefore if 
the s -wave node lies too far from R^ the /' = 2 feature 
becomes larger and narrower, as is seen clearly in Fig. 
8b. A comparison with the experimental width of the 
/' = 2, V = peak leads us to conclude that a consider- 
able flf-wave contribution is present and thus the i'-wave 
node cannot lie to far from R^. This reasoning, however, 
does not tell us on which side of R^ the s -wave node is 
situated. 

We can use the spectra of the higher vibrational levels 
to further constrain the position of the i'-wave. The ratio 
of the purely of -wave /' = 4 peak to the i'-wave compo- 
nent of the /' = 2 peak is proportional to the square of 
the ratio of the ground state wavefunctions at a charac- 
teristic distance R, [30]. A simple estimate of the inten- 
sity ratio of the i'-wave and of -wave contributions to the 



spectral lines can be made based on the approximate 
wavefunctions for the s- and of -waves and is given by: 



d 

s 



k"'R' 



k'Rt 



r'\R, -a)/ (R, - af 



(9) 



where we use the A: ^ expression of J2(kRy)/yk for 
the flf-wave Sindjo(k(Ry — a))l\/k for the i'-wave, and a 
is the scattering length. We can conveniently take Ry to 
be the outer turning point of the OJ v level. An improve- 
ment of the model of the peak ratios involves replacing 
a with the position of the last s -wave node. This follows 
from the modification of the i'-wave wavefunction due to 
the long range — CJR'' potential and is discussed in 
Appendix A. The k dependence shows that, as expected, 
the /' = 4 peaks will disappear for lower temperatures. 

The /' = 2 peak is the dominant feature in the experi- 
mental spectra of the v s 12 vibrational levels. The 
outer turning points of these levels are between 70 Co 
and 200 Co- By Eq. (9) an s -wave node at these internu- 
clear separations would imply a much stronger /' = 4 
peak relative to the /' = 2 peak than observed. We thus 
conclude that there is no i'-wave node between 70 Co and 
200 flo- Since we have already shown that a node too far 
away from R^ leads to an unacceptably small of -wave 
contribution to the v = spectrum, we can also immedi- 
ately rule out a node larger than 200 Oq. Furthermore, a 
small value for the location of the node is also unaccept- 
able as it leads to a of -wave feature that is unacceptably 
weak and a v = 0, /' = 2 level that is unacceptable nar- 
row. Numerical calculations of the peak ratio as a func- 
tion of the shape of the potentials confirm these simple 
arguments. 

Plotting the /' = 2 to /' = 4 peak ratio as a function of 
the shape of the potentials gives the band labeled "peak 
ratio" in Fig. 7. The shape of the potentials at which the 
two bands intersect is the optimal form. Fig. 9 compares 
the theoretical spectra calculated using the best ground 
state potentials with the experiment. The only adjustable 
parameters are the overall height, which is adjusted to fit 
the observed /' = 2 peak and the absolute frequency 
which is adjusted by ~ 2 MHz. The relative peak posi- 
tions and heights are determined from the theory. 

From our final potentials we find zo = 60 Co ± 3 Oo, 
Zi = 73 flo ± 3 flo for the positions of the last s- and 
p-wave nodes, respectively and ai^..., = 52 ao±5 Oq. 
Quoted uncertainties are one estimated standard devia- 
tion (combined standard uncertainty). Other scattering 
properties can be evaluated as well. For example, the 
scattering length 02,2 of two atoms with/a = 2, m/= 2, is 
85 flo ± 3 Oq. This is the scattering length relevant in 
experiments aiming at Bose condensation in doubly po- 
larized samples of Na atoms. 
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Fig. 9. Comparison of theoretical and experimental rotational spectra 
for V = and v = 1 . The theory is scaled to agree with the experimental 
/' = 2 peak height and shifted slightly (~ 2 MHz) in frequency. 

The Na Oi-i scattering length has been discussed in 
the literature before. An experimental measurement of 
fli-i = 92 flo ± 25 flo [31] was based on the thermaliza- 
tion time of a sample with a temperature of 200 (xK. A 
theoretical treatment based on improving on the semi- 
classical RKR potentials with an inverted perturbation 
approach obtained 86-f3 Oo (Ref. [26]). These values are 
consistently larger than our value, although in agreement 
within two sigma if the uncertainties are taken to be one 
sigma. Even without our detailed numerical calcula- 
tions, the observed spectra show that the last/= 2, s- 
wave node cannot lie between 70 Oq and 200 Oq. 



5. Conclusion 

An analysis of the rotational lineshapes in photoasso- 
ciation spectra of the purely long-range Na2 Og' state, 
particularly the lowest vibrational level, places con- 
straints on the possible positions of nodes in the 
3^S(/= l) + 3^S(/'= 1) scattering wavefunctions. By 
combining this information with the known spec- 
troscopy of the Na2 ground states we generate a set of 



potentials which produce scattering phase shifts consis- 
tent with our observed spectra. From the potentials we 
calculate the .?-wave scattering lengths needed as input 
for theories describing Bose condensates. 

Our results reported here are preliminary in that they 
are based on a small data set which limits our ability to 
quantify the effects of the ionizing laser. In future exper- 
iments we plan to acquire a larger data set and also 
investigate spectra in which one or both of the colliding 
atoms are in the 3^S(/'= 2) state. We predict that these 
spectra will be dramatically different from the ones 
reported here and their observation will provide an im- 
portant cross check on the potentials we have derived. 



6. Appendix A. From Nodes to a 
Scattering Length 

This Appendix aims to give an intuitive understanding 
of why for Na2 the last/= 1 , p -wave node zi of the zero 
energy wavefunction lies outside the corresponding 
node of the/= 2 or/= 0, .?-wave. We also relate the 
.?-wave node to a scattering length. 

If we ignore the hyperfine contribution in the multi- 
channel /= 1 and/= 2 Hamiltonians the sole difference 
between the two Hamiltonians is the centrifugal barrier 
£{£+ \)/2iJiR^ where €=\ or 0, respectively. Decreas- 
ing /from one to zero in a continuous fashion makes the 
interaction slightly more attractive and, hence, increases 
the phase that the zero energy wavefunction accumu- 
lates when integrating from R = 0, where the wavefunc- 
tion is zero, to the position of the last /j-wave node. 
Therefore, an .?-wave node lies just inside zi- However 
this does not prove that it is the last .?-wave node. The 
wavefunction could accumulate enough phase in the 
larger R region that it has one more node, i.e., the 
.?-wave potential could have one more bound level than 
the;? -wave potential. The Na hyperfine interaction adds 
small corrections to this picture. This nodal pattern is 
confirmed by full multi-channel close coupling calcula- 
tions for a variety of realistic X'2g and a^'Xu potentials. 
For heavier alkali-metal atoms, however, such a conclu- 
sion need not apply as the hyperfine interaction is larger 
and the centrifugal barrier much lower. We will assume 
that zo stands for the node in ^^\R). The node of the 
/= .?-wave wavefunction, ^J'(7?), is closely related to 

Zo- 

From Sec. 3.3 we know that the wavefunction for 
R>30 ao is described in terms of a single potential and 
the exact wavefunction has a node in this region. This 
allows us to ignore multi-channel complications. The 
connection between the last node and the scattering 
length for the collision between two /a = 1 , rrtf^ = — 1 
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atoms is, therefore, far more tractable. If the atom-atom 
interaction were zero beyond the position of this zero the 
connection is trivial with a\-i = zo- The attractive long- 
range dispersion interaction however is still important. 
In first order the correction to the scattering length due 
to the van der Waals interaction has the form [17,32] 



fli.-i - Zo = Mm-^Yn 






£ 



ARsm{k{R -Zo)) 



phin(^(«-z„)) 






2jU,C( 



{R-z,f 2iJiC6 1 
^^ 30z^ 



(10) 
(11) 



<0. 



For example, if we take the view that Zo = Zi ~ 70 Co Eq. 
(10) implies ai-i = Zo ~ 6 = 64 Oq. A more elaborate 
theory is constructed starting from a zero-energy scat- 
tering wavefunction ij/n{R) as an asymptotic expansion 
in 1/R . The first terms in this expansion can be shown 
to be 



^,lX/^) = (/^-«,,_0.fci-j^.|^ 



,2/A 2 J. 

^^ ^'^ R 



+ 0[(^Cefi, 



(12) 



where ai-i is the scattering length. This wavefunction 
must be zero at zo leading to 



fll.-l = Zo 1 /s^^^n 4^ with Ce = 2lJbC(J^^. (13) 



5 -wave node is at Zo = 37.6 Aq the scattering length is 
infinite or, alternatively, a bound state at threshold has 
appeared. For zo smaller than this critical value another 
node much further out appears. In Fig. 10 the scattering 
length as defined in Eq. (14) as a function of the last 
node in the zero-energy wavefunction is shown. For Zo 
around 70 Co to 80 Co the effects of the — Ce/R "^ are small 
and fli-i = Zo- Near zo = 45 the scattering length be- 
comes negative and for 37.6 Co will become infinitely 
large. 
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Fig. 10. The scattering length versus the last node zo in the zero-en- 
ergy scattering wavefunction. The figure shows the scattering length 
for two assumptions regarding the long-range behaviour of the poten- 
tial. The dashed line corresponds to a zero potential for R> zo and the 
full line corresponds to a — Cf/R^ potential for R>zo- The two 
parameters in the model are the Ce coefficient and the atomic Na 
mass. 
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tra bound 
pure !/«" 
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Vr/-i/4(x 
scattering 
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expression it follows that for zo ~ 42 ao the 
length goes to infinity, or equivalently an ex- 
level appears. According to Ref. [32] for a 
potential the exact ij/uiR) is known analyti- 
a linear combination of yrJinix) and 
) with X = \/2iJicJ¥/{2R^) which leads to a 
length in terms of a zero of the wavefunction 



VC^/-y4(Xo)r(3/4) 

2 Ju4(xo) r(5/4) 



with Xo : 



V2jU,C6/^" 

2zl 



(14) 



where /„(x) is the Bessel function. The scattering length 
as defined in Eq. (14) again has poles, i.e., goes to ± oo, 
as a function of the position of a node in the zero-energy 
s -wave wavefunction. These poles occur at the zeros of 
the function /i/4(xo). For Na this implies that if the last 



The long-range potential is not a pure l/R "^ potential. 
The Cg and higher order terms in the polarization inter- 
action must be included. However, they are small for 
internuclear separations larger than 30 ao- In fact, the 
size of the corrections fall inside the 5 % uncertainty of 
Ce quoted by Ref. [25] and from Eq. (1 1) it follows that 
this adds at most 1 ao to 2 ao to the final uncertainty in 
the scattering length. 
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